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Quantum criticality of metal-insulator transitions in correlated electron systems is shown to belong 
to an unconventional universality class with violation of Ginzburg-Landau- Wilson (GLW) scheme 
formulated for symmetry breaking transitions. This unconventionality arises from an emergent 
character of the quantum critical point, which appears at the marginal point between the Ising-type 
symmetry breaking at nonzero temperatures and the topological transition of the Fermi surface 
at zero temperature. We show that Hartree-Fock approximations of an extended Hubbard model 
on square lattices are capable of such metal-insulator transitions with unusual criticality under a 
preexisting symmetry breaking. The obtained universality is consistent with the scaling theory 
formulated for Mott transition and with a number of numerical results beyond the mean-field level, 
implying that the preexisting symmetry breaking is not necessarily required for the emergence of 
this unconventional universality. Examinations of fluctuation effects indicate that the obtained 
critical exponents remain essentially exact beyond the mean-field level. It further clarifies the 
whole structure of singularities by a unified treatment of the band-width-control and filling-control 
transitions. Detailed analyses on the criticality, containing diverging carrier density fluctuations 
around the marginal quantum critical point, are presented from microscopic calculations and reveal 
the nature as quantum critical "opalescence" . The mechanism of emerging marginal quantum critical 
point is ascribed to a positive feedback and interplay between the preexisting gap formation present 
even in metals and kinetic energy gain (loss) of the metallic carrier. Analyses on crossovers between 
GLW type at nonzero temperature and topological type at zero temperature show that the critical 
exponents observed in (V,Cr)203 and K-ET-type organic conductor provide us with evidences for 
the existence of the present marginal quantum criticality. 

PACS numbers: 71.10.Fd, 71.30.+h 



I. INTRODUCTION 

Metal insulator (MI) transitions in correlated electron 
systems have been a challenging subject of studies for 
decades In this paper, we show from microscopic 
analyses that MI transitions of correlated electrons bear 
an unconventional feature in view of quantum phase tran- 
sitions. This is shown by taking an example of the 
Hartree-Fock approximation of the Hubbard model in 
two dimensions. One might suspect whether the Hartree- 
Fock approximations could yield a phase transition which 
is beyond the conventional scheme. However, we show 
that the simple mean-field approximation allows us to 
get insight into this issue by capturing the correct in- 
terplay of symmetry-breaking and topological characters 
of the transition. We also discuss that the unconven- 
tionality survives beyond the limitation of the mean-field 
study, where the scaling theory [2, 3] formulated for the 
Mott transition is applicable. To make the motivation of 
the present study clearer, we describe somewhat detailed 
introduction to clarify why the present study has a sig- 
nificance and how the unconventionality arises in view of 
quantum phase transitions in general. Then in the latter 
part of the introduction, we review the previous studies 
on MI transitions to position the starting point. 
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A. General perspective in view of quantum phase 
transitions 

A ubiquitous example of phase transitions seen in na- 
ture is the transition between gas and liquid. At the 
critical point of the gas-liquid phase transition, the com- 
pressibility 

1 dV 1 dn 1 dn , . 

~VdP ~ ndP ~ ~rf~d^ ^' 

and the equivalent density fluctuations J dr(n(r)n(0)} di- 
verge, where V and fi are the volume and chemical poten- 
tial, respectively. The opalescence observed in light scat- 
tering is caused by such diverging density fluctuations. 
Nearly a century ago, the critical opalescence found in 
the scale of the light wave length was analyzed by Orn- 
stein, Zernike Q , Einstein and Smoluchovski [f| Q as one 
of the evidences for the existence of atoms. 

The gas-liquid transition is mapped to the ferromag- 
netic transition in the Ising model with the Zi symmetry 
expressed by the Hamiltonian 

H = -J^SiSj-hJ2Sh (2) 

for the local spin Si taking ±1 at the site i, and inter- 
acting with short-ranged neighbors at j and under the 
symmetry breaking field h. The Ising model represents 
the simplest example that phase transitions in general 
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take place because of a spontaneous breaking of symme- 
tries originally hold in the Hamiltonian at h = 0. Al- 
though the gas-liquid transition at a first glance does not 
look breaking the symmetry as in the ordered phase of 
the Ising model, the transition is actually mapped to the 
Z2 symmetry breaking along the first-order coexistence 
boundary. In fact, the up and down spin phases in the 
Ising model are mapped to the gas and liquid phases, re- 
spectively, from which the both transitions turn out to 
belong to the same universality class. In more general, 
thermal phase transitions are all characterized by a cer- 
tain type of spontaneous symmetry breakings. 

Phase transitions by the conventional spontaneous 
symmetry breaking is correctly described by the 
Ginzburg-Landau- Wilson (GLW) scheme UMM- 

In this 

scheme, the free energy F is expressed by a functional of 
a spatially dependent order parameter m(r) after elimi- 
nating other degrees of freedom. F is assumed to allow 
a regular expansion with respect to m(r) such as 



f Oj F^ b 

F = dr[—hm(r) H — -m{r) 2 H -(Vm(r)) 2 H — -m{r) A 

J 2. 2. 4. 

+ •••], (3) 

where b and K are positive constants and the coefficient 
a depends on the control parameter g around the critical 
point g c as a = ao(g — g c )- 

This classical picture of the phase transitions is mod- 
ified in quantum systems. For instance, when the trans- 
verse field — r J2i Sf is additionally applied to the Ising 
model ^ by considering spin- 1/2 quantum spins Sj = 
(Sf,Sf,Sf), T c is suppressed to zero at a critical value 
r c and the transition disappears as illustrated in Fig. [T] 
This may be regarded as a destruction of the symmetry 
broken order by quantum fluctuations instead of thermal 
fluctuations. 




FIG. 1: (color online) Schematic phase diagram of Ising model 
under longitudinal and transverse magnetic fields, h and T, 
respectively. Up and down spin phases are separated by the 
first-order boundary illustrated by shaded(blue) sheet below 
the critical line shown by bold (brown) curve in parameter 
space of T, h and the transverse Ising field I\ Beyond the 
quantum critical point at r c , the transition disappears and 
the up and down spin phases are adiabatically connected even 
through the route on the T = plane. 

In fact, a quantum system may be expressed in the 
path integral formalism, where an additional dimension 



representing the imaginary time or inverse-temperature 
axis is introduced in addition to the d-dimensional spatial 
degrees of freedom. Although the added imaginary time 
direction could be qualitatively different from the spa- 
tial dimensions, one can show that it simply corresponds 
to z additional spatial dimensions, where z is called the 
dynamical exponent. The transverse Ising model is sim- 
ply represented by z — 1. Namely, the quantum phase 
transition at T = in d spatial dimensions realized by 
increasing T to T c may be mapped to the transition at 
nonzero temperature in d + 1 dimensional classical sys- 
tems [l(J. 

However, in more general, if the quantum dynamics 
has a coupling to dissipative gapless excitations such as 
particle-hole excitations in metals, z may become larger 
than unity such as two or three [ll|, [II, [HI- Ferro- 
magnetic and antiferromagnetic phase transitions in met- 
als have been extensively studied from this viewpoint. 
The quantum critical point appears when the critical 
temperature is lowered to zero. Around the quantum 
critical point, low-energy excitations arising from grow- 
ing magnetic correlations have been proposed to cause 
non- Fermi liquid properties (Til [T2I Il3j |. A particular 
type of metamagnetic transitions has also been studied 
around the "quantum critical end point" [3], where an 
apparent symmetry does not change at the critical point, 
whereas an implicit symmetry breaking does occur as in 
the case of the gas-liquid transition. Even in these cases 
with z larger than unity, however, the GLW scheme re- 
mains essentially valid when we consider a proper d + z- 
dimensional system. 

Here the dynamical exponent z is defined from the ra- 
tio of the temporal correlation length £ t to the spatial 
one £ of the order parameter m(r, i) in d dimensions. In 
the case of conventional quantum phase transitions, £j 
and £ are defined by 
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lim 



lim 

t — voo 



log(m(r,0)m(0,0)> 



log(m(0,t)m(0,0)) 



(4) 
(5) 



respectively. More concretely, the dynamical exponent 
^ is defined in the asymptotic limit to the critical point 



9 -> 9c by 



z = lim 



log 6 
log£ 



(6) 



As we will show in the next subsection, in the case of the 
metal-insulator transitions, these definitions should be 
modified, because the correlations in metals decay not 
exponentially but with power laws as a function of dis- 
tance in general. 

In contrast to the description of the phase transition 
by the symmetry breaking, several classes of quantum 
phase transitions at T — are not captured by such ex- 
tensions of the GLW scheme and are not characterized by 
the spontaneous symmetry breaking. One of the simplest 
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examples is found in the transition between a metal and 
a band insulator. The ground-state energy E of nonin- 
teracting electrons near the parabolic band bottom with 
the single-particle dispersion e = Ak 2 in d dimensions is 
obviously described by 

£«^+ 2 )/ d oc^ +2 )/ 2 , (7) 

in the metallic phase for the electron density n and the 
chemical potential fj, measured from the band bottom. 
On the other hand, throughout the insulating phase with 
the chemical potential varied in the gap, E = holds be- 
cause of n = 0. This trivially nonanalytic form of the 
energy as a function of /i clearly represents a quantum 
phase transition with the singularity at fi = 0, but it is 
equally obvious that it does not follow the GLW scheme. 
This transition of noninteracting fermions is not accom- 
panied by any kind of spontaneous symmetry breaking. 
It is rather characterized only by the presence or absence 
of the Fermi surface, which is a topological difference. 
Wen [l5[ has proposed a category of "quantum orders" 
for this class of transitions. 

In this paper, we show that an involved interplay of 
the topological nature at T = with the conventional 
GLW character at T ^ emerges for metal-insulator 
(MI) transitions of correlated electrons. By controlling a 
parameter for quantum fluctuations, a zero-temperature 
critical line characterized by the topological transition 
of Fermi surface from a metal to an insulator switches 
over to a finite-temperature critical line characterized by 
GLW scheme. The phase diagram is now seriously mod- 
ified from Fig. [1] to a new type depicted in Fig. [5] Here, 
the quantum critical point is turned into the starting 
point (or end point) of the quantum critical line which 
represents topological nature of transitions at T = 0. We 
call this starting and marginal point shown as the cross 
in Fig. [2 the marginal quantum critical point (MQCP), 
which is completely different from the conventional quan- 
tum critical point. Then electron correlation effects gen- 
erate an unconventional universality class of phase tran- 
sitions, which is expected neither in GLW category nor 
in simple topological transitions. 

B. Introduction for metal-insulator transitions 

Now we summarize understanding of MI transitions 
for strongly correlated electrons achieved in preceding 
studies in the literature to make clear our motivation of 
the present study. This issue has attracted long-standing 
interest and various important aspects have been clari- 
fied although complete understanding has not been 
reached. MI transitions may occur with coupling to lat- 
tice distortions or be accompanied by effects of random- 
ness in real materials. However, purely electronic origin 
of MI transitions has been a central issue in the strong 
correlation regime, where we find many examples called 
Mott transition. For this purpose, Hubbard model, which 
only considers the local Coulomb interaction on a lattice, 
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FIG. 2: (color online) Proposed phase diagram of MI transi- 
tion for strongly correlated electron systems. A crucial differ- 
ence from Fig. [T]is that the quantum critical line starts along 
the 6 axis (the thick (green) line along a = T = ) beyond the 
endpoint of the first-order boundary called marginal quantum 
critical point depicted by the cross. This difference is a di- 
rect consequence of the topological nature of MI transition at 
T = 0. Here the first-order boundary is again shown by the 
shaded (blue) surface. 

is one of the simplest models to describe MI transitions. 
Hubbard model is defined by the Hamiltonian as 

h = - uAcr c i° + u n ^ n a - m N i, ( 8 ) 

i,j,<r i i 

where cj CT (ci CT ) creates (annihilates) an electron with spin 
a at site i, respectively. The number operator is and 
\i represents the chemical potential. Hopping integral 
between site i and j is represented by iy and U denotes 
the onsite Coulomb repulsion. 

Approximate and simple theoretical descriptions of the 
Mott transition at T = without any symmetry break- 
ing being involved was proposed in seminal works of 
Hubbard [l6| and Brinkman and Rice [l7|. These two 
works offer completely different pictures for the transi- 
tion. Hubbard approximation offers a picture of simple 
splitting of a band into an upper and a lower band with 
a gap between these two bands. The transition to Mott 
insulator occurs simply because the chemical potential 
moves into the gap. The "quasiparticle" weight Z at the 
Fermi level is retained nonzero even at the transition, 
although the "large Fermi surface" with the Luttinger 
volume is not preserved and a small pocket of the Fermi 
surface shrinks to zero at the transition. Brinkman-Rice 
picture obtained from the Gutzwiller approximation, in 
contrast, preserves the Luttinger volume and the Mott 
insulator is realized by the vanishing quasiparticle weight 
Z. Dynamical mean-field theory(DMFT) [n| proposed 
some unification of these two pictures with the forma- 
tion of the upper and lower Hubbard bands reproduced 
together with the vanishing Z in the coherent additional 
"band" . In DMFT results, the Mott transition itself is 
triggered by the vanishing coherent band and the critical- 
ity is ultimately determined from that of the Brinkman- 
Rice scenario. 

Recently, DMFT has been improved by taking into ac- 
count the momentum dependence of the self-energy in 
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the correlator projection method [l9j]. The Mott transi- 
tion obtained from this improvement shows that a transi- 
tion qualitatively different from the original DMFT sce- 
nario occurs, where the renormalization factor Z is kept 
nonzero until the transition [20(. Instead of vanishing 
Z, the transition is realized by shrinkage and vanishing 
of the Fermi-surface pocket with some similarity to the 
Hubbard picture. A crucial difference from the Hubbard 
picture is that the upper and lower Hubbard bands may 
have qualitatively different dispersions accompanie d by 
electron differentiation in momentum space [19, [U, l22j | . 
The vanishing pocket looks to show very anisotropic spec- 
tral weight with strong intensity on the inner-half circle 
and weak intensity on the outer-half circle, making an ap- 
parent arc-type structure for the model on the square lat- 
tice. This arc-type structure is consistent with the exper- 
imental observations in the copper oxides obtained from 
angle resolved photoemission spectra ( ARPES) [23|, [24| . 
This Fermi-surface pocket appears after the Lifshitz tran- 
sition within the metallic phase, where the topology of 
the Fermi surface changes from that of a large Fermi 
surface to pockets. Such Lifshitz transitions seem to in- 
evitably occur, which invalidates the Luttinger theorem 
after the transition to the pockets [2(|. A similar con- 
clusion was obtained by cellular DMFT study for sys- 
tems with a one-dimensional anisotropy [25[ as well as in 
2D (2?| , and by dynamical cluster approximation for 2D 
Hubbard model [261 ] . where the momentum dependence 
of the self-energy was considered in different ways. The 
serious modification of DMFT and the complete change 
in the nature of the MI transition with vanishing pocket 
of the Fermi surface seems to have a deep connection to 
the satisfaction of the hyperscaling and a momentum- 
space differentiation as we discuss below. 

A picture of MI transitions completely different from 
various mean-field approximations including DMFT has 
been proposed from the scaling theory of Mott transi- 
tion [2J, |3[ , where the singularity of the free energy F at 
the Mott transition is assumed to follow the hyperscaling 
form 

Foe £-(<*+*), (9) 

Here, £ in the metallic side is the correlation length de- 
fined by A _1 / d for the metallic carrier concentration X 
and spatial dimensionality d, while it is the localization 
length of carriers in the insulating side. More precise def- 
inition of £ should be given from the correlation length 
of the order-parameter correlation function, which will 
be discussed in the end of this section and will be clari- 
fied that the above relation is correct. The hyperscaling 
is normally satisfied below the upper critical dimension 
d Ul whereas the mean-field theories are satisfied above 
d u . Therefore, the scaling theory and various mean-field 
theories are incompatible in general. 

Before proceeding to more involved discussions on the 
scaling theory, we briefly summarize basic points of crit- 
ical phenomena and universality classes particularly for 
metal-insulator transitions used in the later discussions. 



The universality class of phase transitions is character- 
ized by critical exponents. We first discuss relations for 
critical exponents which is applicable irrespective of the 
validity of cither the hyperscaling or mean-field approx- 
imations. When we take the carrier concentration X as 
the order parameter in the metallic side, the critical ex- 
ponent 5 is defined by 

AcxC 1/5 , (10) 

which measures the scaling behavior of the growth of the 
order parameter as a function of the field £, which is 
conjugate to the order parameter X. Namely, £ is ei- 
ther the chemical potential fi for the filling-control tran- 
sition or U/t in the Hubbard model for the band-width 
control transition. In the latter case, £ is, for example, 
controlled by pressure P in actual experiments. The car- 
rier concentration X is the doping concentration itself for 
the filling-control case while it is the unbound doublon 
(doubly occupied sites) and holon (empty site) concen- 
trations measured from the values at the critical point 
for the band-width control case. 

The exponent v is defined from the relation of the con- 
trol parameter g to £ as 

Z<x\g-9c\- v - (11) 

Here, the control parameter g may also be controlled by 
U/t or /1, which makes some complexity and could cause a 
confusion, because these two quantities can control both 
of C and g. This is somehow in contrast with the simple 
case of Eq. ©, where h and a may be controlled inde- 
pendently, say, by the magnetic field and temperature, 
respectively. In the later microscopic description in Sec. 
IIII1 this somewhat confusing situation will be resolved. 

The exponent 7 is defined from the order-parameter 
susceptibility \ — dX/d( from the scaling 

xoc \g -g c \-\ (12) 

The jump of X at the first-order boundary grows from 
the critical point as 

X^\g~g c f, (13) 

The exponent a is obtained from the second derivative 
of the free energy as 

d 2 F/dg 2 cx\g-g c \- a . (14) 

It should be noted that a for quantum phase transitions 
does not express the exponent for the specific heat as in 
thermal transitions. The set of these exponents has a sig- 
nificant meaning because the set specifies the fundamen- 
tal nature of the phase transitions called the universality 
class. 

Conventional scaling laws such as Widom's law 

00-1) = 7 (15) 
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and Rushbrooke's law 



a + 2(3 + 7 = 2 



(16) 



may be satisfied for these exponents. So far, the control 
parameter g is not fixed yet and one may have several 
choices. Some of the exponents need trivial transforma- 
tion depending on the definition of g 28] . 

When the hyperscaling assumption is employed, the 
length scale which diverges toward the critical point is 
unique, and it must be the mean carrier distance given 
by 



(17) 



Then the free energy expressed by Eq. @ leads to 

F oc X {d+Z ^ d . (18) 

This is the form at £ = and g— g c — 0. Scaling relations 
are derived by further adding the term of the "chemical 
potential" £ and the term of coupling to g — g c as 



F = -(X + B Q {g - g e )X+ + CX {d+z ^ d , 



(19) 



where Bo and C are constants and 4> should be sepa- 
rately given from a physical consideration or microscopic 
derivation of the transition. Note that the definition of g 
is now fixed as the quantity to control the transition at 
C = 0. The scaling between X and £, namely, Eq. (JTUJ) is 
obtained by putting g — g c = in Eq. (|19|) and minimize 
the free energy by dF/dX = 0. This yields 



5 = z/d. 



(20) 



The scalings between X and g, namely, Eqs. (fl"2"|) and 
(fTB")) are obtained by putting £ = in Eq. (JTHJ) and 
minimize the free energy by dF/dX = 0. This yields 



and 



z — d 



1 z+d(i-4>y 

From (fTTj) . (fT5|) and (|17p. a further relation 

v = f3/d 



(21) 



(22) 



(23) 



is obtained. Note that the scaling laws Eqs. (| 1 5[) and 
(fT6]) are obviously satisfied for Eqs. JSJ-fSH]). From the 
scalings ©, CP and ([17]), we obtain F oc \g - g c \ u[d+z) 
for nonzero g— g c . From this and Eq. (114[) . the Josephson 
relation 



a = 2 - v(d + z) 



(24) 



is obtained. In addition, Drude weight is shown to fol- 
low Q 



To estimate accurate exponents of microscopic mod- 
els, in principle, one needs to go beyond the mean-field 
approximation, because the reliability of the mean-field 
approximation itself has to be critically tested. Avail- 
able theoretical tools for this purpose is severely limited, 
because of various difficulties in theoretical approaches 
for strongly correlated electrons. In the literature, there 
exist only small number of available estimates from un- 
biased approaches such as the quantum Monte Carlo 
and path-integral renormalization group methods. For 
the filling-control transition of Hubbard model or the 
t-J model on a square lattice, by using the quantum 
Monte Carlo method [H, [H HU and power Lanczos 
method (32j, the charge compressibility has been shown 
to follow the scaling dX/dfi oc A" 1 , from which 8 is es- 
timated from Eq. (jTOJ) as 5 ~ 2. The result obtained 
from the path-integral renormalization group method is 
also consistent with this scaling [35]. When the hyper- 
scaling holds, 8 = 2 and Eq. (j2T))) in two dimensions 
leads to z = 4. On the other hand, simple GLW scheme 
always gives S > 3 and is incompatible with the numer- 
ical results. In fact, the exponents of the Ising model 
are 5 ~ 4.8 for three-dimensional systems and S = 15 
in two dimensions. Any other types of symmetry break- 
ings indicate S > 3. Note also that the degeneracy tem- 
perature (effective Fermi temperature) Tp is scaled by 
T F oc oc X z l d resulting in T F oc X 2 for z = 4 and 
d = 2, if the hyperscaling is satisfied 0]. 

Experimentally, the exponent 5 has been examined 
for filling-control transitions by measuring the chemical 
potential shift with increasing doping by photoemission 
studies in several high-T c copper oxides (39l . |40| . In par- 
ticular, I^S^CaC^Os+j, and (La,Sr)2Cu04 show pin- 
nings of the chemical potential at low doping, indicating 
an enhancement of the charge compressibility and con- 
sistency with (5 = 2. 3 He adsorbed on a graphite surface, 
offering a unique two-dimensional fermion system with 
strong correlation effects, appears to show an anomalous 
suppression of Tp near the commensurate solid phase in 
accordance with the trend of Tp oc X 2 consistently with 
z = 4 41] with the hyperscaling being satisfied. 

In addition, the Drude w eig ht has been calculated to 
be consistent with D oc X 2 [H, [13] which suggests z = 4 
from the above relation (j2l)|) if the hyperscaling is satis- 
fied. These two independent estimates of z coincide each 
other and support the hyperscaling. 

Other exponents are directly related to g dependence 
of F and it is fixed when <f> is specified. If one put z = 4 
and d = 2, 
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2(3-0)' 
1 



a = 2-3- 



3-. 



D oc X 1+i - z - 2 ^ d . 



(25) 



ft 



1 



(26) 
(27) 
(28) 
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and 
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(29) 



where we have used relations for the equal-time structure 
factor S and susceptibility x, given by 



are obtained. 

This unconventional universality class has further been 
analyzed and (f> = 2 is obtained from a plausible assump- 
tion on the quasiparticle dispersion [3y, H3, HI] . Then, 
the universality class is proposed to be characterized by 
a set of the exponents 

z = 4,a= -1,P= 1,7 = 1,6 = 2,u = 1/2 and 77 = 0. 

(30) 

At finite temperatures, it has been theoretically pro- 
posed that the MI transition is equivalent to the gas- 
liquid transition and its universality class belongs to the 
Ising universality class [12, S3[ ■ Experimentally, from re- 
cent careful study of V2O3 by Limelette et al [4J|, they 
have obtained the critical exponents which are consis- 
tent with that of the 3D Ising model, namely, (3 ~ 0.34, 
5 ~ 5.0 and 7 ~ 1.0 in the region very close to the crit- 
ical point whereas they have obtained consistency with 
the mean field exponents for the Ising model in the most 
part away from the critical point. On the contrary, recent 
experimental results on k-(ET) 2 Cu[N(CN) 2 ]C1 0, S3 
indicate the critical exponents f3 = 1,7 = 1 and 6 = 2, 
which are consistent with Eq. (|30[) and indicates an evi- 
dence for the relevance of the present universality. 

Before closing this section, we remark on the defini- 
tion of £ and the dynamical exponent z. In the case of 
MI transitions, correlation functions do not decay expo- 
nentially as a function of distance in the metallic phase. 
Therefore, the conventional definition of £ given in Eq. 
((4]) is not applied. For example, noninteracting electrons 
have a density-density correlation function (X(r)X(O)) 
which asymptotically decays as 1/r 3 with Friedel oscil- 
lations in two-dimensions [301 ]. Such an asymptotic de- 
cay crossovers to a constant at shorter distance, where 
the crossover length is determined from the inverse of 
Fermi wave number fop. Since 1/r 3 decay is absolutely 
convergent in the integration, the integration over dis- 
tance is contributed dominantly from the region inside 
of the length scale of kp . From this, £ is given by 
£ oc k~ l 



X 1 / d . More precisely, the correlation length 



£ is defined by 



rf/ /dr(A(r,0)A(0,0)) 

5 - V (x(o,o)x(o,o)) • ^ ) 

In the similar way, the temporal correlation length £ t is 
defined by 

Jdt(X(0,t)X(0,0)) 
§ * (X(0,0)X(0,0)) ' { ' 

Then the dynamical exponent z is obtained from 



S = J dr{X(r,0)X(0,0)) af d 
X = [ drdt(X(r,t)X(0 7 0))^^- d 



leading to 



(34) 
(35) 

(36) 



z = hm -j , 

g^gc log £ 



(33) 



C. Outline of the present work 

In this paper, we show that the basic structure of the 
phase diagram for the MI transition involving the quan- 
tum critical line and the boundary for the first-order 
transition is correctly captured by the Hartree-Fock the- 
ory of the Hubbard model. A remarkable result in this 
paper is that the Hartree-Fock approximation is able 
to derive the above unconventional universality class at 
MQCP from a microscopic derivation of the free energy. 
The obtained universality class turns out to be consistent 
with the scaling theory. The mechanism of the emer- 
gence of MQCP and resultant unconventional universal- 
ity is clarified from the present microscopic study. It is 
ascribed to the selfconsistent positive feedback between 
metallic carrier density and the amplitude of an insulat- 
ing gap already developed in metals. More concretely, 
increase of the preformed gap amplitude decreases the 
carrier density to lower the kinetic energy, which in turn 
further enhances the preformed gap. On the contrary, 
the decreasing preformed gap allows more carriers, which 
further reduce the gap to gain the kinetic energy of car- 
riers. These nonlinear effects cause the first-order transi- 
tion and MQCP emerges at the border between the first- 
order transition and the quantum critical line. Further- 
more, the Hartree-Fock theory turns out to remain essen- 
tially correct even when one goes beyond the mean-field 
approximation and takes into account quantum fluctua- 
tions. This is because the transition always occurs as that 
at the upper critical dimension irrespective of the spatial 
dimensionality. We further show that the Hartree-Fock 
solution has a finite-temperature crossover to the conven- 
tional Ising universality class. We discuss implications 
and consequences of this unexplored type of universality 
class emerging near MQCP sandwiched by GLW critical 
and quantum critical lines. We argue that low-energy 
excitations are described neither by the original single- 
particle picture nor by the bosonic density fluctuations 
when the both characters of transitions meet. 

We also discuss comparison of the present result with 
the experimental indications obtained so far. It is shown 
that the critical exponents at MQCP together with the 
finite-temperature crossover are completely consistent 
with the whole aspects of the experimental results. 
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A part of the discussions in this paper is already briefly 
given 49]. In this paper, we present the results of quan- 
tum MI transitions in greater detail. In particular, we 
show details of the procedure to obtain the non-GLW 
free-energy expansion and critical exponents together 
with thorough discussions on implications. 

The organization of this paper is as follows: In Sec. 
im we show details of our Hartree-Fock approximation. 
In Sec. IIII1 we present Hartree-Fock solutions of the 
MI transition by introducing the preexisting order pa- 
rameter. It captures some essence of the unconventional 
universality class of the transition from metals to insula- 
tors caused by the electron correlation effects. Using the 
free-energy expansion, we analytically obtain the critical 
exponents of quantum MI transitions, which belongs to 
an unconventional universality. We also discuss implica- 
tions beyond the mean- field approximations in Sec. IIIII 
In Sec. IIV1 we show numerical results of Hartree-Fock 
approximation at T = 0, which are completely consis- 
tent with the analytical results in Sec lIIIl In Sec. [V] 
finite-temperature crossover is calculated. It shows a 
well-defined crossover from the Ising universality in the 
close vicinity of the critical point to the quantum univer- 
sality governed by MQCP when the parameter deviates 
from the critical point. Section IVT1 describes comparisons 
of the present implications with relevant experimental re- 
sults. Section lVIIl is devoted to summary and discussions. 

II. HARTREE-FOCK APPROXIMATION 

In this paper, to capture the essence of quantum MI 
transitions, we extend Hubbard model ([S]) by considering 
the nearest-neighbor repulsion term proportional to V on 
the N s = L x L square lattice. It is defined by 

ff = ^(£i(k) + &(k))4 CT Cfc CT 

k.CT 

+uj2nitnu + vj2 N i N i-^J2 Ni ' ( 3? ) 

i (ij) 1 

where Ni = (n^ +ri;|) are number operators. We restrict 
the hopping terms to the nearest neighbor pairs given by 
the dispersion £i (k) = — 2i(cos fc^+cos k y ) and that of the 
next-nearest neighbor pairs by £2(k) = At' cos k x cos k y , 
respectively. 

Using this extended Hubbard model, we study quan- 
tum MI transitions in two-dimensional system within the 
Hartree-Fock approximation. Although the Hartree-Fock 
analyses have been done for the Hubbard model in the 
literature [13, 13, the criticality of MI transitions has 
not been examined. By considering the MI transitions in 
the ordered phase, we are able to capture a remarkable 
aspect of MI transitions from a microscopic model. 

To consider the MI transitions within the Hartree-Fock 
approximation, we consider a symmetry broken long- 
range order such as charge or antiferromagnetic (AF) or- 
dering with the periodicity commensurate to the lattice. 



Then we study the MI transitions in the symmetry bro- 
ken phase. We note that essence of the MI transitions 
does not depend on the type of the ordering. For the 
extended Hubbard model defined in Eq. (|37p . two sim- 
ple orders become stable. When AV > U, the extended 
Hubbard model has a tendency to charge ordering, where 
the empty and doubly occupied sites alternatingly align 
within the mean-field level. The order parameter to of 
checker-board-type doublon alignment is defined by 

™ = jr^2i( n n +^i)exp(iQri)), (38) 

S i 

where Q is the ordering wave vector (n,Tt). For U > 
AV ', the antiferromagnetic ordering becomes stable. The 
order parameter of the antiferromagnetic ordering toaf 
is defined by 

m A F = jf- 53 ((rift ~~ ex P(*Q r i))- ( 39 ) 

S i 

Hereafter, we mainly consider the charge ordering. 
However, within the mean-field approximation, the anti- 
ferromagnetic order and charge order play the same role 
on MI transition and by replacing the order parameter 
m with the antiferromagnetic one, one can easily reach 
the same results for the nature of the MI transition. 

Here we take the mean field as 

, , , , (n + mcxpfoQi-j)) 

(n iT ) = ) = - , (40) 

with n being the average charge density given by 
^/N s J2i( n i] + n il)- Using this mean field, we decouple 
the interaction term. 

The on-site interaction term is decoupled as 

i i 

= yEK(« + -' Qri )) 

_ UN s (ni+ m i) 
4 

The nearest-neighbor interaction term is decoupled as 

(ij) (ij) 

= inVY^Ni-AmV^Nie-KpiiQri) 

i i 

- 2VN s (n 2 - to 2 ). (42) 

Finally we obtain a Hartree-Fock Hamiltonian(iJyp') 
in the momentum space as 

H H F = 53 ( £ ( k ) + ~Y + AnV ~ A C ka C k<r 
k,(T ^ ' 
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where we define e(k) = £i(k) + ^(k). 

Diagonalizing the Hamiltonian leads to two bands of 
the form 



S ± =e 2 (k)±^i(k) 2 + A 2 -M, 



(44) 



where A = mg and we now take g = AV—^ as the control 
parameter, because in this section we mainly consider the 
band- width control transition. In Eq. (I44p . we dropped 
a constant term n(g + U). 

Using this Hartree-Fock band dispersion we obtain the 
free energy 

T 



Frf = -jj-tog Z HF 

= -^E lo s( 1 + ex P 

S k,7j V 



E v (k) 
T 



U(n 2 + to 2 ) 



+ 2V{n 2 - to 2 ) , 



(45) 



where the suffix rj takes ±. 

From this free energy, the order parameter m is deter 
mined from the self-consistency condition 

l = ^ v /(E_(k))-/(£ + (k)) 

TV. 



.9 



k 



£L(k) 



(46) 



where f(x) is the Fermi-Dirac distribution function. The 
particle density n is given by 

' ^/(25_(k)) + /(25+(k)). (47) 



N. 



E 

k 



Details of numerical calculations of Eqs. (|4"6"]) and (|4"T|) 
are shown in Appendix |XJ 

It turns out that the self-consistent equation (|4"6")1 with 
Eq. (|4"4"]l for the charge ordering is equivalent to the 
self-consistent equation for the AF ordering defined in 
Eq. ((39)) ji?], Hg|. The equivalence holds by replacing 
the charge order parameter to with the antiferromag- 
netic(AF) order toaf and by putting g — {7/2. However, 
one should be careful about this mapping. This map- 
ping is complete only on the Hartree-Fock level. The 
charge ordering at the (it, it) periodicity is the conse- 
quence of the discrete symmetry breaking while the AF 
ordering is realized by the continuous symmetry break- 
ing of SU (2) symmetry. Therefore, if it would be exactly 
solved in two-dimensional systems, the charge ordering 
can indeed exist at finite temperatures, though AF or- 
dering cannot exist at finite temperatures by Mermin- 
Wargner theorem. Therefore, the Hartree-Fock approxi- 
mation captures some essence of the charge ordering at fi- 
nite temperatures, while the AF ordering at nonzero tem- 
peratures in two dimensions contains an artifact of the 
Hartree-Fock approximation, although the development 
of the charge gap itself is not an artifact, which leads to 
the Mott gap formation even without any long-range or- 
der. The physics of unconventional quantum criticality, 
which we clarify later, emerges from the gap formation 
itself and therefore the essence is captured in the Hartree- 
Fock calculation. 



III. FREE-ENERGY EXPANSION AND 
CRITICAL EXPONENTS 

In this section, we analyze singularities of the free en- 
ergy around the MI transition within the Hartree-Fock 
approximation at T = and obtain the critical exponents 
of the transition. When the order parameter to and the 
gap A = mg are sufficiently small near the charge-order 
or antiferromagnetic transition, the insulating gap does 
not open for nonzero positive t' . This is because, the 
upper band E + and lower band _E_ overlap for a small 
gap. From Eq. [|4"Ij). the top [bottom] of the lower [upper] 
band is located at (k x ,k y ) = (tt/2,tt/2) [(71", 0)], respec- 
tively. The continuous MI transition occurs when the top 
of the lower band reaches the same energy as that of the 
bottom of the upper band as shown in Fig. [3] Therefore, 
the condition of the continuous MI transition is given by 



£+(7r,0,A c ) = J B_( 7 r/2,7r/2,A c ) ! 



(48) 



where A c is the critical gap of the continuous MI transi- 
tion. From this, at T = 0, we obtain the critical gap A c 

as 



A r = 2t'. 



(49) 



Metal 



Transition point 




FIG. 3: (color online) Quasiparticle band structures of the 
Hartree-Fock solution for metallic phase (left) and the con- 
tinuous MI transition point (right) are shown. Broken (blue) 
line shows the free-particle band dispersion defined by e(k) = 
— 2i(cos k x + cos k y ) + 4t' cos k x cos k y . The definition of the 
upper (lower) band is defined in Eq. (|44[) . Fermi energy (if is 
represented by the horizontal line. In the metallic phase, the 
bottom of the upper band at M point ((71-, 0) and its equivalent 
points) is lower than the top of the lower band at (tt/2, tt/2) 
and its equivalent points so that the indirect insulating gap 
is closed. 

This quantum MI transition occurs between charge- 
ordered metal and charge-ordered insulator, thus the 
symmetry does not change at the transition. It is le- 
gitimate to expand the free energy around nonzero order 
parameter measured from the critical point of the MI 
transition. One might speculate that the free energy ex- 
pressed only by the order parameter A in Eq. (|45[) would 
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be regular around the MI transition. However, as we will 
see below, this expansion has a piecewise analytic form in 
the metallic and insulating phases separately with a jump 
of the expansion coefficient at the transition point. This 
somewhat unconventional form of expansion is derived 
from the fact that some of the expansion coefficients ex- 
plicitly depend on the density of states at the Fermi level, 
which has a discontinuity at the MI transition. This free- 
energy structure with a jump is in a complete contrast 
with the GLW form. 

In this section, we consider the band-width control 
transition at half filling (n = 1) in the canonical ensem- 
ble. In the scheme of the canonical ensemble, we add 
a term /in to the free energy defined in Eq. (|45|) . The 
Hartree-Fock free energy Fqhf is given by 



(50) 

where rj runs over + and — . Hereafter, we describe Fchf 
as F. 



A. Free-energy expansion 

Since A dependence of F is explicitly given in Eq. ([50]) . 
it allows a formal expansion of F in terms of A around 
a nonzero value A c , by taking an expansion parameter 
Y = A c - A. The region Y > (Y < 0) represents the 
metallic phase (insulating phase) , respectively. The form 
of expansion is given as 



F(A) = F(A C )+AY+^Y 



AY* 



where the coefficients are given by 



OF 











A=A C 





(51) 
(52) 

(53) 



B 



d 2 F \ 

5A 2 ! 

1 v 



N > k 



A=A C 



0/(£„(k,A)) 



A 



dA 



w 2 



A,. 



c = 



d 3 Fi 

~dA3\A=A c 

1 ^( d 2 f(E n (k,A)) 

N s 



dA 2 



with the definitions 



W 1 =v 1 , 

W 2 =v x -A 2 c v 3 , 

W 3 = v 3 A c - v 5 A 3 c . 



(54) 

3VF 3 , 
(55) 

(56) 
(57) 
(58) 



Here, v n is defined by 

w ™ = -TT E^i 2 + A^/OMk, A c )). 



(59) 



k . r} 



To derive Eqs. ([54")) and (p)5"|) , we neglect terms those are 
proportional to 



a/(^(k,A)) d a 



dA 9A^ 



(60) 



A,. 



because of the following reason: At T = 0, this term has 
nonzero value only at £"77 (k, A) — (on the Fermi sur- 
face) because df(E ri {k,A))/dA becomes the delta func- 
tion. At the MI transition point, since the Fermi energy 
has the same value as that of the top [bottom] of the lower 
[upper] band, df(E_(k, A))/dA [df(E + (k, A))/flA] 
has nonzero value only at (tt/2, 7r/2)[(7r, 0)], respec- 
tively. However, at the MI transition point, namely 
at (tt/2, tt/2) and (77, 0), £1 = is satisfied, then 

d/dA ^A/ A 2 ^ becomes zero. From this, the term 

defined in Eq. ([617]) vanishes at the MI transition point 
A = A c at T = 0. 

Next, we derive more explicit forms of coefficients B 
and C at T = 0. We first consider the contribution of 
the first term of Eq. fl5I]) and (USD at T = 0. The first 
term of Eq. (jB"4")) is given by 



0/(£„(k,A)) A 



— V 



dA 



J dED v (Er,)G(T,En 




(61) 



where D V (E V ) is the density of states of the upper (lower) 
band for 77 = +1(— 1), respectively, A is a cut-off and 
G(T, E) is defined by 



G(T,E) 



eME/T) 



T(l + exp(E/T)) 2 ' 



(62) 



We introduce a constant a for the linear coefficient be- 
tween A c — A and /1 as 



= <5(A C - A) - 2t'. 



(63) 



Detailed derivation of a at T = is described in Ap- 
pendix |5] Explicit form of a at T = is given by 



2t' - \/t 2 - V 



2t' + yjt 2 - t' 2 



(64) 



The function G(T, E) becomes the delta function at T = 
0. Then, Eq. (|6T|) is discontinuous at the MI transition at 
T = 0, because the density of states at the Fermi level is 
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discontinuous between the metallic and insulating sides. 
In fact, at T = 0, we obtain the relation 



= - (D+{1 + a) + D-(l -a)) for Y > 
= for Y < 



where D + (D_) is the density of states at the bottom 
(top) of the upper (lower) band, respectively. Explicit 
forms of D_|_ and Z?_ are given by 

1 



Da 



D- 



4xf 



1 



2n\/t 2 - t' 2 



(65) 
(66) 



Detailed derivation of D + and _D_ are given in Appendix 

m 

Using the similar derivation, we obtain the representa- 
tion of the first term of Eq. (|5"5"|) as 



y ^ dED rj (E ri 



dG(T, E v ) 



dA 



A,-. 



From this, we obtain the relation as 



(67) 



= -(1 + a) 2 R+{T) + (1 - a) 2 i?_(T) 
+ %^(l + 5)(3 + a) 
+%^(l-a)(3-a) forF>0 
= for y < 

where R±(T) and Q±(T) are defined by 

i?„(T)= [ A dEll ^^±G(T,E v ) (68) 



Q„(T) = / dE, 

J-A 



D v( E v) 



-G(T,E n ). (69) 



Finally, at T — 0, we obtain the free-energy expansion 
for metallic and insulating sides as 



F = AY + ^Y 2 + ^Y 3 . 



F = AY HlY 2 + ( 4Y 3 . 



Bi_ 

21' '3! 



where 



A = -(--W 1 )A c 
9 



(70) 
(71) 

(72) 



and the coefficients of the metallic side (A < A c ) are 
given as 

B m = -W 2 - ((1 + &)£>„ + (1 - a)Di) + - (73) 

9 

C m = -3W 3 - R+(0) (1 + &f + R-(Q) (1 - a) 2 

+ 9A {1 + fi)(3 + &) + RM. (i _ a)(3 _ & ), 

(74) 



and the coefficients of the insulating side (A > A c ) are 
given as 



B, = -W 2 



Ci = -3VFs 



1 



(75) 
(76) 



We emphasize that, while A is common, the coefficients 
B and C in the region Y > have values different from 
that in the region Y < at T = 0. Although the coef- 
ficients B and C jump at Y = 0, F is regular (namely, 
piecewise analytic) within each region Y > and Y < 0. 
Then the Taylor series expansion (|70| as well as (17T1) is 
allowed in each sector Y > and y < separately and 
they are unique. The coefficient Ci is negative if t' is 
not too large, because the sign of W3 is determined from 
the dominant contribution proportional to V3 , which has 
the positive sign due to the contribution from the lower 
band dominating over the upper-band contribution. On 
the other hand, C m is positive for the relevant region of 
t', namely for not too large value of \t'\, because the dom- 
inant contribution is then from the term proportional to 
Q+. Therefore, the MI transition takes place in the pa- 
rameter regions of C m > and C, < 0, which makes F 
always bounded from below within the expansion up to 
the cubic order of Y. In fact, we will show below that 
the expansions up to the cubic order is enough for the 
purpose of deriving critical properties. Note that this 
expansion up to the cubic order is consistent with that 
obtained in the scaling theory and its phenomenological 
extension for two-dimensional systems [36|, H3, HI] • We 
will show below that the universality class obtained in 
the present paper by the microscopic Hartree-Fock cal- 
culation indeed agrees with the prediction of the scaling 
theory and further reveals a microscopic mechanism of 
the emergence of the unconventional universality. 

The carrier density X , which is the sum of the electron 
density A_ and the hole density X + in the metallic phase 
near MI transitions, is given as 



X = X 4 



X. 



1 



7rA, 



(i + 5)y. 



(77) 



Details of the derivation of Eq. ([77)1 are given in Ap- 
pendix A. Since X is proportional to Y in the metallic 
phase, the singularity of X is the same as that of Y . 
Therefore, whichever of X or Y is taken as the order 
parameter, the critical exponents are the same in the 
metallic side. 



B. Critical exponents 

Using the free-energy expansion obtained in the previ- 
ous subsection, we are able to obtain the critical expo- 
nents of the quantum MI transitions. 

Phase diagram of the MI transitions in the A-B m plane 
is shown in Fig. @]for the region close to A = B m = 0. 
Metallic states are realized when the minimum is at a 
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First-order line 

F 



Q\ Y 
\ 
i 
\ 





l m 5=2 


Metallic Phase 




'F 


\ 




i 


1 



FIG. 4: (color online) Phase diagram in A-B m plane. The 
first-order transition line is given by the condition A ~ 
SB m /8C m with B m < 0, which separates the metallic phase 
A < 3B m /8C m from the insulating phase in the other side. 
For B m > the transition becomes continuous and the tran- 
sition line is given by A = (see Eq. ((Ml)). Typical Y 
dependence of the free energy along the first-order transition 
line, in the metallic and insulating phases are also shown in 
the insets. The definition of the free energy is given in Eqs. 
(|70|) and (|71[) . The real Y dependence of the free energy is 
represented by the solid lines, where the metallic (insulating) 
side corresponds to the region Y > (Y < 0), respectively. 
Broken lines represent the behaviors of Eqs. (|7U|I and (J7TJ ex- 
trapolated to the other sides Y < and Y > 0, respectively. 
For the derivation of the phase boundary, see the text. 



positive Y m , whereas insulating states are stabilized when 
the minimum is at a negative Yi. From Eq. (|51|) . a local 
minimum appears at 



Y — 

1 m. — 



B rn + Va» 2AC n 



(78) 



if B 2 > 2AC m and at 



Yi = 



-B % + ^B 2 - 2AC l 



(79) 



if Bf > 2ACi is satisfied. Because F is bounded from 
below, there always exists the absolute minimum either in 
the region Yi < or Y m > 0. The location of the absolute 
minimum and the phase boundary shown in Fig. [5] are 
clarified below by classifying the parameter values into 
several cases. Then the critical exponents are derived in 
each region. We remind that B t > B m and Cj < < C m 
are always satisfied in all the relevant regions. 



1. Quantum critical line 

The continuous MI transition line (quantum critical 
line) is determined from the condition A = and B rn > 0. 
Along this quantum critical line, the free-energy mini- 
mum exists at Y = Q. For nonzero A near the quantum 
critical line A = 0, Eqs. (|75|) and (|T9")) are reduced to 



and 



Y„ 



Y 



A 

B m 



A_ 

~Bi 



(80) 



(81) 



respectively. Because the coefficient A plays the role of 
the conjugate field to Y, the critical exponent 5 is defined 
from the growth of Y for small but nonzero A as 



Y-A 1 ' 6 . 



(82) 



Therefore, S = 1 is obtained in both of the metallic and 
insulating sides. 

The susceptibility in the metallic and insulating phases 
are given by 



X 



' d 2 F y 1 


1 


AY) 2 ) 


B m 




1 


\d{Yf) 


" B~i 



(83) 



respectively. Along the quantum critical line, B is of 
course always positive for the region B m > in whichever 
of the insulating and metallic sides. Therefore, the sus- 
ceptibility does not diverge and the critical exponent de- 
fined by 



is given as 



7 = 0. 



(84) 



(85) 



As we mentioned in Eq. (fH)) in Sec. QJ the critical 
exponent a does not represent the singularity of the spe- 
cific heat in the case of the quantum phase transition. 
It rather characterizes the singularity of \ ~ d 2 F/dA 2 , 
because in the present case around the quantum critical 
line, A oc g — g c holds. Near the quantum critical line of 
the MI transition, the free energy F at the minimum is 
scaled by A as 



F oc A 2 , 



(86) 



as is shown by substituting Eqs. (|80|) and (JSTJ) into Eqs. 
([TO")) and CI]). From Eq. dHJ) , this indicates a = 0. The 
exponent (3 is not well defined on the quantum critical 
line [5l| . because it is not possible to control the param- 
eter g to cause the phase transition without applying the 
field A conjugate to the order parameter. 
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The free energy form (|70f and (|7Tj) with nonzero A, Bi 
and _B, m is equivalent to the case of the transition between 
metals and band insulators for noninteracting electron 
systems, if one replaces Y with the carrier density A 
of noninteracting electrons. By utilizing the equivalence 
with this noninteracting system, in the present case of 
the quantum critical line, we also obtain £ oc A -1 / 2 oc 
\g — g c \~ x ^ 2 leading to v — 1/2. On the other hand, one 
can also calculate the dynamical susceptibility \ of the 
carrier density in the noninteracting system, which yields 
x/S oc Ep 1 oc kp 2 oc £ 2 . Here, E F is the Fermi energy. 
This leads to the dynamical exponent z = 2. 

Now we summarize the critical exponents on the quan- 
tum critical line: a = 0,7 = 0,(5 = l,v = 1/2 and z = 2. 
This is common to both of the insulating and metallic 
sides of the transition. 



2. Marginal quantum critical point 

When B m becomes zero, the quantum critical line ter- 
minates at the marginal quantum critical point (MQCP) 
and the first-order transition boundary starts. Around 
MQCP, we show that the exponents depend on the ap- 
proach to the critical point either from the metallic side 
or from the insulating side. To clarify such route depen- 
dence, we introduce four routes to approach MQCP as is 
illustrated in Fig. [5] 



Insulator 




(d) Bn 



FIG. 5: (color online) Four routes to approach MQCP. Solid 
line in the negative B m region represents the first-order transi- 
tion line. The quantum critical line is given by the condition 
A = and B m > 0(see text). Route (a) is the way to ap- 
proach MQCP by controlling "the external field" A from the 
metallic region. The singularity of the Y along the route (a) 
is characterized by the critical exponent 8. Route (b) is the 
way to approach MQCP from the insulating phase. Route (c) 
is the way to approach MQCP along the first-order transition 
line A ~ §B, 2 n as is derived in Eq. flMJ. Route (d) [(d')] is 
the way to approach MQCP along the quantum critical line 
from the metallic [insulating] phase. 

From the free-energy expansion, the coefficient A rep- 
resents the conjugate field to the order parameter and B 
represents the control parameter g. Therefore, we obtain 
the critical exponent S through the route (a) or (b) in the 
following: 



Route (a) 

In the metallic side, through the route (a) in Fig. [5j the 
coefficient A approaches zero by keeping B m — 0. From 
Eq. (|75|) . neglecting B m , we obtain the singularity of Y 
as 



Y 

- 1 ID 



'2\A\ 



oc 



On 



cc\A\i 



(87) 



Therefore, the critical exponent 6 in the metallic side is 



6 = 2. 



(88) 



Note that this exponent is consistent with the value ob- 
tained in the scaling theory 0, H, [H, H3, 1381 as well as 
with that in the quantum Monte Carlo study [29|, [3(1 HH . 

The dynamical exponent of the present transition can 
be calculated from the correlation functions of the order 
parameter X. Namely, from the equal-time structure fac- 
tor S for the carrier density X defined in Eq. (|34|) and 
the susceptibility x defined in Eq. (|35|) . z is given from 
Eq. ©. By calculating S and x, we obtain z = 4 [501 ] . 
Alternatively, we can derive z from the following argu- 
ments: Characteristic length scale of the MI transitions 
is given by the mean carrier distance £, which diverges 
at the continuous MI transition point including MQCP. 
Then, in two dimensions, carrier density X is propor- 
tional to £ -2 . The equal-time structure factor S has the 
same scaling oc £ 2 by its definition. On the other hand, 
from <5 = 2, the susceptibility of the carrier density x 
is scaled by A -1 oc £~ 2 . Then the dynamical critical 
exponent z is determined from the condition 



(89) 



as z — 4. It also yields that £ has the singularity £ oc 
\g — g c \ v with v = 1/2. The large dynamical exponent 
2 = 4 is consistent with the previous predictions obtained 
beyond the mean- field approximations 0, H, [3(1 [3l| . 

The origin of such an unusually large z is now clar- 
ified from the present study: The dynamical exponent 
normally expresses the dispersion of quantum dynam- 
ics. Namely, one might expect that the carrier dynamics 
would be described by the dispersion e oc k 4 . However, 
this contradicts the quasiparticle dispersion (|4"4"|) , which 
is generically proportional to k 2 . This puzzle is solved 
when we realize that the gap amplitude decreases with 
the increase of carrier number. Because of this reduction 
of the gap, when the carrier density is fixed, the upper 
(lower) band edge as well as the Fermi level becomes 
lower (higher) than that we might expect from the rigid 
band picture obtained from the quasiparticle dispersion 
of insulators. Because of this selfconsistent reduction of 
the gap, the Fermi level does not move as quickly as k 2 
with increasing A and is rather pinned near the origi- 
nal band edge even when the carrier number increases. 
This results in an effectively flat dispersion and actually 
a higher exponent with z = 4 comes out. This is an intu- 
itive picture for the filling-control transition. For the case 
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of the band-width control, where the chemical potential 
is fixed rather than the filling, the reduction of the gap 
further allows more electrons and hole carriers, which is 
again interpreted by an effectively flatter dispersion. 
Route (b), (d') 

In the insulating phase, through the route (b) and (d') 
in Fig. only A becomes zero while Bi remains nonzero 
and positive even at MQCP. This condition generates 
the same criticality as that of the quantum critical line. 
Therefore, we obtain the critical exponents a = 0, 5 = 1, 
7 = and z — 2. 

Route (c) 

Next we consider the critical exponent (3 defined by 



Y ~ IB I' 3 

1 jump l-L'ml , 



where Yj ump is the jump of Y across the first-order tran- 
sition line located in the region B m < 0. Near MQCP, 
from Eqs. (|78j) and (f79|) . Yj ump is given by 



Y = Y 

1 jump — 1 r 



B m \X A 

Cm Bi 



(91) 



where we define A = 1 + y/B^ - 2AC m /\B m \. From this 
definition, it turns out that A is of order unity. On the 
first-order transition line, the free energy in the metallic 
phase has the same value as that of the insulating phase, 
namely 



F{Y m ) = F(Yi). 
Equation (|92]) leads to 

C m (X\B m \f 



A\\B m \ B m (\\B m \) 2 



2C 2 



(92) 



A 2 

= -2M- 



Near MQCP, \C m \ and \Bi\ are always nonzero and have 
amplitudes of the order unity. We assume that, for small 
quantities A and B m , A has the order e\ and B m has the 
order From Eq. (|93|) . the relation between ei and e 2 
is given by 

0(e 1 )0(e 2 ) + 0(e|)~0(e?). (94) 
, the relation between ei 



From this, to satisfy Eq 
and £2 should be 



(95) 



Substituting Eq. into Eq. ||7S} leads to A = 3/2. 

Namely, on the first-order transition line near MQCP, 
the relation between A and B m should be 



A 



8G n 



(96) 



Using this relation, near MQCP, the singularity of Yj ump 
is given by 



Y 



■SBl 



&C m Bi 



3 1 -B n 

2Cri 



oc \B n 



(97) 



Therefore, the critical exponent /? is given by 

13= I. 



(98) 



Route(d) 

The route (d) along the quantum critical line (^4 = 
and B m > 0) always keeps A = 0. Although the sus- 
ceptibility x remains finite on the quantum critical line 
itself (x — 1/B m ), the amplitude l/B m increases toward 
MQCP and it diverges in the limit to MQCP because of 
B m — > 0. The singularity of the susceptibility 



(90) is given by 



X = IS - 9c\~ 



d 2 F 



d{Yf 



(99) 



(100) 



\B m \- x 



« \g-g c \ 



Therefore, in the metallic phase, the susceptibility di- 
verges at MQCP with the critical exponent 



7 = 1. 



(101) 



Near MQCP, the free energy F at the minimum in the 
metallic side is scaled by the control parameter B m as 



FtxBt 



(102) 



For MQCP, B m oc g — g c holds in contrast to A oc g — g c 
for the quantum critical line. Then from Eq. (|14[) . we 
obtain the exponent a = — 1. We remark again that a 
has no direct relation with the singularity of the specific 
heat. In the insulating side, as mentioned above for the 
route (d'), a = is obtained. 

Now the critical exponents for MQCP are summarized 



a = -1,(3 = 1,7 = 1,5 = 2, is = 1/2 and z = 4 (103) 
for the metallic side and 

a = 0, 13 = 1,7 = 0,5= 1,1/ = 1/2 and z = 2 (104) 



for the insulating side. The scaling laws (| 1 5|) and IT6l) 
as well as the hyperscaling Josephson relation (|2~4"|) asso- 
ciated with the hyperscaling Eq. © are all satisfied in 
both of the two sides. 



C. Crossover of the critical exponent 5 

In the metallic region, Y grows from zero at MQCP 
when A and/or B m deviate from zero. As a function of 
the distance from MQCP defined by r = \J A 2 + B 2 t in 
the A-B m plane, Y grows as Y oc r for A = because 
of Eq. (|9"T|) whereas it is scaled as Y oc for B m = 
because of Eq. {37]) • Then a crossover between Y oc r 
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and Y <x r 3 must exist. Here, we specify the location of 
this crossover boundary below. 
Near MQCP, Y m is given as 



Y 

- 1 if 



-B m + y/B* - 2AC„ 



— —B' 



BL 



AL 



(105) 
(106) 



where we define B' m = B m /C m and A' m = — 2A/C m > 0. 
In the region where the condition A' m 3> B' m 2 is satisfied, 
Ym behaves as 



Y 

- 1 it, 



(107) 



This indicates that MQCP exponent(<5 = 2) dominates in 
this region. In contrast, in the region where the condition 
AL <C B' 2 is satisfied, Y m behaves as 



Y m ~ -B' m + B' m (l 



BL 



A' 
_ 

B 



r«l4 



(108) 
(109) 



From this, in this region the criticality of the quantum 
critical line represented by S = 1 dominates. In Fig. [6l 
the crossover line is represented by the broken line. 



zero temperature. Figure [7] shows numerical estimate of 
B m near MQCP from Eq. ([75|. 

From this, we estimate the location of MQCP at 
t'Jt = 0.05571, g c = 0.616180 and t'Jt = 0.36455, gjt = 
1.38144, respectively. We call MQCP at t' c = 0.05571, 
MQCPi and MQCP at t' c = 0.36455, MQCP 2 , respec- 
tively. 
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FIG. 7: (color online) Coefficient B m on the line A = 
0(A = A c ) near MQCP. We estimate the points B m = 
as t' /t = 0.05571, and 0.36455. In the positive B m region, 
the continuous MI transition occurs, while in the negative 
B m region, the first-order MI transition occurs. 



E. Finite temperature effect 
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FIG. 6: (color online) Crossover line between the scaling re- 
gion Y oc r and V oc r', which is represented by the broken 
(blue) line defined by A' ~ \B' m \ 2 . 



D. Location of marginal quantum critical point 

From the free-energy form, the location of MQCP is 
determined by the condition 



A = and B m = 0. 
The line A = is determined by 
A = A,, 



(110) 



(111) 



with the definition A c = 2t'. From Eq. (jllip . the rela- 
tion between g and t' for A = is given, because A is an 
implicit function of g and t'. We calculate B m from Eq. 
(175)1 with the condition Eq. (TTTTj) . We then find that B m 
becomes zero at two points, namely two MQCPs exist at 



In this subsection, we analyze how the free energy 
and critical properties are modified at nonzero temper- 
atures. Since the discontinuity of the coefficient arising 
from the singularity of the density of states is immedi- 
ately smeared out by the Fermi distribution at T > 0, 
the free-energy expansion should convert into the expan- 
sion of the conventional tfi 4 theory in the critical region. 
The universality class of MQCP then switches over to the 
Ising universality class. Here, we show how the expan- 
sion in Eq. I|5ip breaks down and the Ising critical region 
appears. 

The discontinuity of the coefficient B and C comes 
from the integration of the form 

,A 

P= / dEJ2D TI (E)H(E)G(T,E), (112) 

where H(E) is a slowly varying function of E in the range 
E ~ 0(T). Now the density of states D + (E) is nonzero 
only in the interval 

- (1 - oc)Y < E < A, (113) 

whereas D_(E) is nonzero in the window 

-A<£< (l + a)Y. (114) 

The correction factor proportional to a arises from the 
chemical potential shift from the case of Y = 0. By 
considering Eq. (f6"2"]) , Eq. (|112p is reduced to 



P = 2H(0) A ; (l + tanh 



(1 - r)a)Y 
2T 



), (H5) 
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if the energy dependence of the density of states D v (E) 
can be ignored. Since (1— a)/2 and (l+a)/2 in Eq. qil5[) 
are the quantities of the order unity, from Eq. (|112p , for 



Y ^> T, P behaves similarly to that of the step func- 
tion as shown in Fig. [8j Therefore, for Y ^> T, the 
expansion Eq. ([5Tj) is valid and the MQCP criticality 
appears. We call this region the quantum region. On the 
other hand, for F<T,P behaves differently from that 
of the step function as shown in Fig. [8j and behaves as 
a smooth function of Y. Therefore, the expansion (|5"Tj) 
breaks down and the conventional GLW expansion Eq. © 
becomes justified, where the Ising criticality appears. We 
call this region the classical region. We expect that fluc- 
tuations modify the mean-field exponents for the Ising 
universality because the upper critical dimension of the 
Ising universality is four. We do not go into detail on 
this point. 

Quantum region ~Y»T xP 

vr 




Ising region ~y«r 

Metal y 



FIG. 8: (color online) Schematic behavior of P as function 
of Y . In K > T region P behaves as the step function and 
the expansion (|51[1 is valid, while in F < T region P behaves 
as Fermi distribution function and the expansion (|51[) breaks 
down. 

Now, we clarify width of the quantum region for each 
route shown in Fig. [9] 




Metal ^ t' 



FIG. 9: (color online) Routes to the critical point in plane of g 
and t' at fixed temperature T. Solid thick (red) line represents 
the first-order transition line and it terminates at the finite- 
temperature critical point represented by the (yellow) circle. 
Route (c) is the way to the finite-temperature critical point 
along the first-order transition line with constant T. Routes 
(a) and (b) are the routes to the critical point with fixed tem- 
perature from the metallic and insulating phases, respectively. 
Broken line represents the crossover line, where \dA/dg\ has 
the maximum value. 

Route (a) 



In the quantum region, the singularity of Y m is given 
from Eq. ([57]) as 



Ym oc 




a\9~ 9c\\ 



where a is a (/-independent constant given by 



a = lim 



IS— 9o|->0 \g — g c \ 2 ]j C, 




(116) 



(117) 



From this, the condition of quantum region Y 3> T is 
jiven by 



, T 2 
9~9c\ > — ■ 



(118) 



by 



Route (b) 

In the quantum region, the singularity of Yi is given 



Y oc — - b\g - g c \, 



where & is a ^-independent constant defined by 

1 A 



lim 



Is-ScHo \g - g c \ Bi 



(119) 



(120) 



From this, the condition of quantum region Y 3> T is 
given by 



T 

\9-9c\ > j- 



(121) 



Route (c) 

In the quantum region, the singularity of Yj ump is given 
from Eq. (J57J as 



v \3> 

* jump f t 



cit'-a 



where a ^-independent constant c is defined by 

1 B m 

c = lim 



|t'-t;,H0 \t' - t' c \ C„ 



(122) 



(123) 



From this, the condition of quantum region 7 > T is 
given by 



T 

c 



(124) 



IV. NUMERICAL HARTREE-FOCK RESULTS 

AT T — 

In this section, we show numerical results of MQCP 
obtained by numerically solving the Hartree-Fock self- 
consistent equation (|46|) itself without relying on the free- 
energy expansion (|51|) . We will show that numerical re- 
sults are consistent with the free-energy expansion in the 
previous section. 



16 



A. Phase diagram at zero temperature 

We first present the phase diagram at T — in Fig. [TO] 
We find two MQCPs in the phase diagram in agreement 
with the free-energy expansion. 
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FIG. 10: (color online) Phase diagram at T = 0. Broken 
(solid) line represents the first-order (continuous) transition 
line, respectively. Near MQCP, charge- ordered metal appears. 
Detailed phase diagram near MQCP is shown in Figs. Illf a) 
andfTTfh) 



At t' = 0, the perfect nesting occurs and the charge 
ordered phase begins from g — 0. For a narrow window 
below the triple point (t' ~ 0.085i, g ~ 0.694i), charge- 
ordered metal(COM) appears between paramagnetic- 
metal(PM) and charge-ordered insulator(COI) phases as 
is shown in Fig. Illf a) in detail. The triple point repre- 
sents the point where PM, COM and COI coexist. The 
order of the transitions from PM to COM is continuous, 
while the order of transition from COM to COI changes 
from the continuous to the first order with increasing 
t' /t. The continuous transition line between PM and 
COM terminates at the triple point. The quantum crit- 
ical line between COI and COM terminates at MQCPi 
(f ~ 0.05571*, g ~ 0.6161*). 

For large t' region (*' > 0.35*), another triple point 
appears at *' ~ 0.357i, g ~ 1.366£ and COM appears 
between COI and PM as is shown in Fig. [TTT b) in de- 
tail. The quantum critical line between COM and COI 
terminates at MQCP 2 (£' ~ 0.3645i, g - 1.381£) again. 

Figure [T2] shows the jump of Y(Yj ump ) as a function of 
t' across the first-order MI transition line at T = 0. The 
(red) squares and (purple) triangles show Yj ump between 
COM and COI and between COI and PM, respectively, 
as a function of t'. Yj ump vanishes at MQCP. 
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FIG. 11: (color online) (a) Phase diagram near MQCPi (t' ~ 
0.05571t, g ~ 0.61614). The continuous (first-order) transi- 
tion is shown by the solid (broken) line, respectively. MQCP 
is shown by the circle, where the first-order transition line 
terminates. The cross point represents the triple point where 
COI, COM and PM coexist, (b) Phase diagram near MQCP 2 
(t' ~ 0.3645t, g ~ 1.38W). 
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FIG. 12: (color online) Yj ump as a function of t' on the first- 
order MI transition line. For the first-order transition between 
PM and COI shown by triangles, Yj U m P is slightly larger than 
A c = It' . 



B. Critical exponents of MQCP 

We clarify the singularity of Y near MQCP numeri- 
cally. Hereafter, we only consider MQCP at t' c — 0.36455, 
namely MQCP2 through the route (a)~(d) in Fig. [5] 
Critical exponents of MQCPi is the same as that of 
MQCP2- We show that the numerical Hartree-Fock re- 
sults are consistent with the results obtained by the free- 
energy expansion in the previous section. 

Route (a) 

First, we confirm that the critical exponent 5 in the 
metallic region defined by 



Y m cc \A\ s and Y„, 



a\9 - 9c 



(125) 



is consistent with the value obtained by the free-energy 
expansion. We perform least square fitting for Y m as 
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a function of g c — g near MQCP. Figure [13] shows the 
results of the least square fitting. We estimate the critical 
exponent 8 as 



8 = 2.0001(6). 



(126) 



This is consistent with 8 = 2 obtained by the free-energy 
expansion. We estimate the coefficient a defined in Eq. 
CUD as 



13.48(4). 



(127) 
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FIG. 13: (color online) Log-log plot of Y m as function of g c — g 
at t' /t = 0.36455. Solid (green) line shows the result of least 
square fitting by assuming the function y = ax 1 ^ 6 , where 
y — Y m and x = (g c — g). Results of least square fitting are 
given as 1/8 = 0.49996(2) and a = 13.48(4). Critical gap A c 
is obtained as A c = 2t' = 0.72911t and the critical interaction 
g c is estimated as g c = 1.38144. 

Route (b) 

We also obtain the critical exponent 8 in the insulating 
region. The definition of 8 is given by 

Y z (x\A\^ and Y m ~ b\g - g c \i (128) 
The numerically obtained critical exponent 8 is 

5 = 1.000000(4). (129) 
The coefficient b defined in Eq. (|120|) is estimated as 

6 = 0.996429(4). (130) 

Now it was confirmed that the Hartree-Fock result is con- 
sistent with the exponent 8=1 obtained by the free- 
energy expansion. 
Route (c) 

Next, we consider the critical exponent (3 defined by 



Yj um p OC l-Sml oc c\t t c \ 

We estimate the critical exponent (3 as 
(3 = 1.01(1). 



'1/3 



(131) 



(132) 



This critical exponent is consistent with f3 = 1 obtained 
by the free-energy expansion. The coefficient c defined 
in Eq. Q123[) is estimated as 



c= 35.2(1). 



(133) 
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FIG. 14: (color online) Crossover points of the critical expo- 
nent 8 determined from g — g c dependence of Y for t' = 0.3650 
near MQCP2 at t' = 0.364555. Away from the transition 
point, behavior of Y for t' = 0.3650 asymptotically ap- 
proaches that of Y for t' = 0.364555. Crossover point is de- 
termined from the crossing point of a\g — gc] 1 ^ 2 and a'\g — g c \. 
The t' dependence of the crossover point is shown in Fig. 1151 
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FIG. 15: (color online) Crossover boundary between 8 = 2 
region and 8 = 1 region in plane of g — g c and t'. The inset 
shows log- log plot of crossover in plane of |p — gr c | and t' c —t' ■ 
Solid (green) line shows the result of least square fitting near 
MQCP by assuming the function y = kx 2 , where y = \g — g c \ 
and x = \t' c — t'\. 



C. Crossover of the critical exponent 8 

The critical exponent (5 has different values between 
the quantum critical line and MQCP and it may show a 



18 



crossover when the critical regions meet. Here, we show 
numerical results for the crossover of the critical exponent 
6 at T = 0. As already mentioned, at MQCP, singularity 
of Y is well fit by the function a\g — gd 1 ^ with 5 = 2 
and a = 13.48. Away from MQCP, 5 crosses over from 2 
to 1 near the quantum critical line. In Fig. I14( we show 
the least square fitting for g c — g dependence of Y by 
assuming the function a'\g — g c \ x ^ ', where 5 = 1. The 
crossover point is estimated by the crossing point of the 
function a\g — gd 1 ^ 2 and a'\g — g c \, which is illustrated in 
Fig. EH The inset in Fig. [TBI shows that the t' dependence 
of the crossover point is well fit by the function \t' — 1' 2 \ 2 . 
This is consistent with the analytical result. 



V. NUMERICAL HARTREE-FOCK RESULTS 

ATT/0 

To understand quantitatively, we numerically study 
how critical exponents of MQCP cross over to those of 
Ising critical points by finite-temperature effects. As 
shown in the previous section, the free-energy expansion 
predicts that crossovers of the critical exponents depend 
on the routes (a)~(c). 

Route (a) 

In the metallic phase, we show the crossover of the 
critical exponent 5 near the critical point for an example 
of T = O.OOlf in Fig. [TBI Here, 6 clearly crosses over from 
2(MQCP) to 3(Ising mean-field value). We estimate the 
width of Ising region by the crossing point of two solid 
(green and blue) asymptotic lines. The value of g at the 
crossing point is denoted as g* . The Ising region is given 
by W < W* for 



W=\g c -g\ (134) 
. In the example in Fig. [111 the 



and W* = \g c - g 
crossover occurs at 



W* 



0.6 x \Q- 7 t. 



(135) 



On the other hand, the free-energy expansion in Eq. 
(|118p indicates that the crossover is given by 



T 2 



0.001 
13.48 



t ~ 10" 8 i, 



(136) 



where a is estimated in Eq. (|127p . They are consistent 
each other. 

For several choices of temperatures, the same analy- 
ses were performed and we determined the Ising region. 
We show the result in Fig. fTTT a) . In the free-energy ex- 
pansion, the crossover W* is expected asymptotically as 
W* ~ T 2 at low temperatures. The scaling of W* as a 
function of T appears to show somewhat smaller expo- 
nent p than two for the fitting W* ~ T p . The origin of 
this discrepancy is likely to come from temperature de- 
pendence in a. Except for this discrepancy, the order of 
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FIG. 16: (color online) Log-log plot of Y m as function of 
\g — g c \ near the critical point at T = O.OOlt. Two thin solid 
(green and blue) lines show the results of least square fitting 
by assuming the function y = ax 1 ^ and y = a'x 1 ^ 2 , respec- 
tively, where y = Y m and x — g — g c . We obtain a = 0.924(1) 
and a' = 14.701(3). 



(a) Metallic phase 
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(b) Insulating phase 
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FIG. 17: (color online) (a) Phase diagram of crossover bound- 
ary between Ising and the quantum regions for S in the metal- 
lic phase plotted in plane of W/t and T/t. (b) Phase diagram 
of crossover boundary between Ising and quantum regions for 
6 in the insulating phase plotted in plane of W/t and T/t. 



the Ising region estimated from the Hartree-Fock calcula- 
tions is consistent with that of the free-energy expansion. 

Route (b) 

In the insulating region, we show the numerically ob- 
tained phase diagram for the crossover in Fig. I17f bh To 
obtain the crossover point from 6 = 3 for the Ising region 
to 6 = 1 for MQCP, we have used a criterion different 
from the metallic side in the route (a). The slope of Y as 
a function of W is not a monotonic function of W and it 
crossovers first from S = 3 to an even smaller value, af- 
ter which it increases to 5 = 1. To specify the crossover 
point, we have taken it as the point which shows the 
minimum slope in the crossover region. For example, at 
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T = O.OOli, W* is estimated as 

W* = 2.23 x 1Q-H. 



1.6x10 



-3 



(137) 



The free-energy expansion indicates from Eq. (|12ip that 
the crossover point W* for 5 is given by 



FE 



T 
~b 



O.OOli 
0.99 



10 _3 i 



(138) 



where b is estimated in Eq. (|130[) . As far as the order of 
W, it indicates that the Hartree-Fock numerical results 
are consistent with that of the free-energy expansion. 
Figure [TTT b) is obtained by performing the same analy- 
ses for several choices of temperatures. The temperature 
dependence of the crossover is estimated as W* ~ T 11 . 
This is rather consistent with W* ~ T obtained from the 
free-energy expansion. It is likely that a slight discrep- 
ancy between two results comes from the temperature 
dependence of b. 
Route (c) 

We also discuss the crossover of the critical exponent 
in the route (c). For example, near the critical point at 
T = O.OOli, crosses over from l(MQCP) to 1/2 (Ising 
mean-field) in a monotonic fashion. The crossover W* is 
estimated by the same method as the route (a) as 



W* = 0.00037- 10~H. 



(139) 



In the free-energy expansion, the crossover is estimated 
from Eq. (fl~24|) as 



W, 



T O.OOli 



FE 



32.4 



10" 4 i, 



(140) 



where c is estimated in Eq. (|133|) . The numerical 
Hartree-fock result is consistent with Wp E . 

Figure [T8l shows W at the crossover for several choices 
of temperatures. The temperature dependence of Ising 
region is estimated as W* — T 7 . which has somewhat 
smaller slope than W* ~ T expected from the free-energy 
expansion. Again, it is likely that the discrepancy comes 
from the temperature dependence of c. Numerical results 
suggest that we need to consider the overall tempera- 
ture dependences of a, b and c on the quantitative level. 
However, we do not go further into details on this point. 
At least the existence of two critical regions governed 
by MQCP and Ising classes is now well established with 
crossovers, which are qualitatively consistent with the 
free energy expansion. The numerical results show how 
the asymptotic scaling for the location of the crossover 
derived in the free energy expansion are modified in re- 
gions distant from the critical point. 



VI. COMPARISON WITH EXPERIMENTAL 
RESULTS 

In real materials, strictly speaking, well defined Mott 
critical point has been observed so far only at finite tem- 
peratures, as we discuss below. There, Ising critical expo- 
nents should be observed in the region sufficiently close to 
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FIG. 18: (color online) Crossover boundary for the critical 
exponent plotted in plane of W/t and T/t. 



the finite-temperature critical point. However, because 
the Ising region is not wide enough at low temperatures 
and crosses over to MQCP criticality, quantum critical 
exponents may be observed in actual experiments. Al- 
though we have to be careful in comparing our results 
for the very simple model with possible artifact of the 
Hartree-Fock approximation with experimental results 
on more complicated systems, it may offer an insight 
into the relevance of the quantum criticality. Here, we 
compare our results with the experimental results of k- 
(ET) 2 Cu[N(CN) 2 ]Cl reported by Kagawa, Miyagawa and 
Kanoda [46| . The critical exponents in the metallic phase 
show a clear evidence of the crossover in our numeri- 
cal results for 6 and 0. At T c = 0.01i(see Fig. flTTa)). 
the classical Ising region is observed for W < 10 _6 i. 
This critical temperature is comparable to the exper- 
imental value T c — 40K with the estimated transfer 
— 0.1- 0.2 eV obtained from the Hueckel calculation [52| . 
In k-(ET) 2 Cu[N(CN)2]C1, l.OMPa roughly corresponds 
to O.OOli [52|, [53], which is estimated from the pres- 
sure dependence together with the dependence on the 
anion species in the experimentally obtained phase dia- 
gram. The smallest distance from the critical point ac- 
cessible in the experimental resolution is estimated as 
\P - P c \ ~ 0.1MPa~ 10~ 4 i, which is much larger than 
10~ 6 i suggested from the theoretical estimate above. 
This indicates that, in the experiments, the distance be- 
tween the true critical point and the available point clos- 
est to it is two orders of magnitude larger than the dis- 
tance of the crossover boundary from the critical point. 
This is consistent with the experimental result that the 
singularity of the conductivity is well fit only by the quan- 
tum critical exponent 5 = 2. Other exponents in this 
k-ET compound shows 7 = 1 and = 1. These are com- 
pletely consistent with the present MQCP universality. 

We have also obtained that 8 in the insulating phase 
shows a crossover. The critical exponent 6 changes into 
S = 3 (Ising) from 5 = 1 (MQCP). We estimate the Ising 
region as W < W* - 10~ 3 i for T c = O.Olf. This W* 
corresponds to about l.OMPa by assuming that l.OMPa 
corresponds to O.OOli. The Ising region in the insu- 
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lating phase is much wider than that of the metallic 
phase. From this, we have more chance for observing 
the Ising region (5 = 3) in the experiments. The exper- 
imental result shows that effective charge gap A closes 
with the singularity as A ~ (P c — p)°- 4 ± 01 i n the in- 
sulating phase. This indicates that 6 = 3 ~ 2. This 
critical exponent may be consistent with the Ising mean- 
field value 5 = 3. If one were able to approach the Ising 
critical region (\P — P c \ -c l.OMPa) in the organic con- 
ductor of k-ET compound, since the system is quasi-two- 
dimensional, the 2D Ising critical exponent 5 — 15 should 
be obtained. It could show a crossover from the quan- 
tum exponent 5 — 1 governed by MQCP to the Ising 
mean-field value S = 3 in the region closer to the critical 
point with an extremely narrow critical region of 5 = 15. 
However, since the experimental data are obtained in the 
parameter space not close enough to the Ising region of 
the critical point, the mean- field critical exponent may 
be obtained similarly to the case of (Vi-zCr^C^ [44J. 
Although the pressure dependence of the parameter value 
for the transfer is so far not available, high critical tem- 
perature ~ 450K of (Vi-zCr^Os likely enables the ob- 
servation of the Ising criticality. 

As we mentioned in SecU some of the high-T c cuprates 
such as P^S^CaC^Og+y and (La,Sr) 2 Cu04 show con- 
sistency with 5 — 2, whereas some of other compounds 
such as (Ca,Na)2Cu02Cl2 shows a quicker shift, rather 
consistently with 5=1. This variety suggests that pa- 
rameter values of various copper oxides correspond to a 
diverse range from the region close to MQCP to that 
deeply in the quantum critical line. The normal state 
properties as well as superconducting transition temper- 
atures may depend on the distance from MQCP, which 
is an intriguing issue left for the future. In the case of 
filling-control transitions, the first-order transition ap- 
pears as an intrinsic instability to the phase separation. 
The phase separation is transformed to the electronic 
inhomogeneity with finite-size domains under the con- 
straint of the long-range Coulomb interaction and the 
charge neutrality. The distance from MQCP directly 
measures the tendency for the instability toward such 
inhomogeneities. The diversity of the cuprate supercon- 
ductors in terms of this tendency can be classified by this 
distance. 

3 He adsorbed on a graphite surface offers a unique 
purely two-dimensional fermion system. In addition to 
the anomalous suppression of Tp that is consistent with 
2 = 4 [4l| as mentioned in SeclH recent detailed study at 
lower temperatures has revealed richer structure in the 
specific heat [lij . The second layer of He approaches the 
registered phase (namely, a solid phase with commen- 
surate periodicity to the first-layer solid at the density 
4/7 in terms of the first layer density) with increasing He 
pressure. Near but below the density at the registered 
phase, the lower temperature structure of the specific 
heat suggests that the large residual entropy expected 
from the large dynamical exponent z — 4 is released be- 
low lOmK. A strong crossover around 10% doping be- 



tween low and high density regions is suggested from the 
transfer of the entropy release, where the release from the 
specific-heat peak structure at 10-100mK for the lower 
density (more than 10% doping) region is replaced with 
the release at lower temperatures around or less than 
lmK for the higher density (less than 10 % doping). The 
release at this low-temperature region may arise from the 
superposition of the contribution from carrier motion and 
spin, which are originated from two different regions of 
Brillouin zone, respectively, as analyzed by the electron 
differentiation in momentum space with a simplified two 
fluid picture (36j . This transfer of the entropy-releasing 
temperature implies a drastic topological change of the 
Fermi surface in the ground state such as Lifshitz transi- 
tion to the small Fermi pocket or approximate arc struc- 
ture. If a clear Lifshitz transition occurs, we expect a 
singular dependence of thermodynamic properties on the 
density measured at low temperatures below lmK. An- 
other conceivable case is that such a singular dependence 
is not observable when the topological change occurs in 
the pseudogap region with faint spectral weight out of 
the " arc" . This may be analyzed by measuring the data 
at temperatures lower than lmK in more detail around 
this crossover, which is an intriguing issue left for future 
studies. If Lifshitz transition occurs, the transition to the 
registered phase occurs also as a topological transition 
with vanishing Fermi pockets to the registered phase. So 
far, it is not clear whether the transition to the registered 
phase is eventually realized across the quantum critical 
line or near MQCP, because the lower temperature data 
below O.lmK is not available. 

Physics of electron differentiation tightly associated 
with the present unusual universality of metal-insulator 
(or liquid-registered phase) transitions is a common cen- 
tral issue of the high-T c cuprates and monolayer 3 He. 
The overall experimental indications show the relevance 
of the unusual quantum criticality around MQCP. More 
detailed analysis is left for future studies. Studies on the 
crossover effect beyond the Hartree-Fock approximation 
are also left for future studies. 



VII. SUMMARY AND DISCUSSION 

Quantum MI transitions between symmetry broken 
metals and insulators have been studied within the 
Hartree-Fock approximation using the extended Hubbard 
model. By increasing the onsite Coulomb interaction U 
and nearest-neighbor Coulomb interaction V, this model 
at half filling has a tendency to charge ordering or anti- 
ferromagnetic ordering within the mean-field level. Pre- 
existing gaps generated by such spontaneous symmetry 
breakings allow us to capture the essence of electron cor- 
relation effects on metal insulator transitions. In this 
paper, we have mainly considered the case of the charge 
ordering in the region 4V > U . However, it is essentially 
the same for the antiferromagnetic ordering in the region 
U > 4V and the present study can be easily extended, 
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where the criticality does not change. 

Using the Hartree-Fock free energy, we have studied 
the criticality of MI transitions. To clarify the nature 
of MI transitions, we have performed the free-energy ex- 
pansion as 

F A = F(A C ) + AY + |r 2 + |y 3 , (141) 

in terms of the gap amplitude measured from the MI 
transition point , Y = A c — A. This free-energy expan- 
sion has an unconventional feature, where the coefficients 
B and C have different values between B m and C m in 
metals and B t and Cj in insulators with jumps at the 
transition point Y = 0. Although the expansion is not 
regular at Y = 0, Y dependence of F is regular in a 
piecewise analytic way in each metallic ( Y > ) and 
insulating ( Y < ) region separately, which makes the 
expansion unique. Because of the jump from C m > 
to Ci < at Y — 0, the whole expansion (|141[) up to 
the cubic order of Y is bounded from below and the free- 
energy minimum appears at a finite Y . Such jumps of the 
coefficients at Y = clearly violate Ginzburg-Landau- 
Wilson scheme of phase transitions and indicate that the 
MI transition at T = is not explained by the concept of 
the spontaneous symmetry breaking. The origin of the 
jumps is ascribed to the jump of the density of states 
from a nonzero value in metals to zero in insulators in 
two dimensions. It reflects the topological nature of the 
MI transition, which is caused by the disappearance (or 
emergence ) of the Fermi surface without a change in the 
symmetry. 

From the free-energy expansion in Eq. p41[) , we find 
that both continuous and first-order MI transitions oc- 
cur at T = 0. The numerical estimates for the relevant 
parameter region of the extended Hubbard model show 
Bi > and B t > B. m , whereas C\ < < C m . When 
B m > 0, the region A < represents the metallic phase, 
while the insulating state is stabilized for A > 0. For 
B m > 0, the MI transitions are always continuous. The 
quantum critical line with the continuous MI transition is 
determined by the condition A = 0. When B m becomes 
zero, the quantum critical line terminates at the marginal 
quantum critical point (MQCP)[v4 = 0, B m — 0] and the 
first-order transition line begins in the region B m < 0. 
Our calculated coefficients A and B m specify the loca- 
tions of MQCP at t' = 0.05571 and t' = 0.36455 for the 
extended Hubbard model. 

In Sec. 4, we have studied the MI transitions by solv- 
ing the self-consistent equation numerically. The ob- 
tained phase diagram in t'-g plane at T = contains 
two MQCPs. We have confirmed that the numerical re- 
sults are consistent with the analytical results obtained 
in Sec. 3. 

Our mean-field exponents at MQCP is essentially exact 
beyond the mean-field theory. Substituting these critical 
exponents to the Ginzburg criterion [9J, we obtain the 



upper critical dimensions d c as 

a c = 2 = 2. 

v 

Since the two-dimensional system is at the upper critical 
dimension, the mean-field critical exponents of MQCP 
should be exact except for possible logarithmic correc- 
tions. 

Using the free-energy expansion in Eq. (|14ip . we ob- 
tain the critical exponents of MQCP as a = — 1, 0=1, 
7 = 1, S = 2, v = 1/2, r\ = and the dynamical expo- 
nent z = 4 in the metallic side of the transition. These 
critical exponents show that MQCP belongs to an uncon- 
ventional universality class different from those obtained 
from phase transitions caused by the spontaneous sym- 
metry breaking. 

The discontinuities of the coefficients B and C ex- 
ist, in the strict sense, only at T = 0. At nonzero 
temperatures, the jump is immediately smeared out, 
which makes a crossover to a different universality. We 
find that the first-order MI transition boundary extends 
to nonzero temperatures and terminates at the criti- 
cal point. Around the finite-temperature critical line, 
the free energy follows the GLW scheme. The crit- 
icality is perfectly categorized by the Ising universal- 
ity. This means that the universality class of the finite- 
temperature transition may be characterized as a symme- 
try breaking transition in contrast to T = transitions. 

The peculiarity of the universality at MQCP is as- 
cribed to the fact that MQCP appears as the marginal 
point between the conventional GLW transition at 
nonzero temperature with emergence of the spontaneous 
symmetry breaking, and the topological transition at 
T = 0, each of which are described by already known 
physics, while its connecting point is not. This unconven- 
tional universality class of MQCP is consistent with re- 
cent experimental results of k-(ET) 2 Cu[N(CN) 2 ]C1. Al- 
though we have focused on the band-width control MI 
transitions, these critical exponents are the same as those 
of the filling-control MI transitions. 

We have clarified how the crossover between Ising 
and MQCP universalities appears. It has been shown 
that even the finite-temperature critical point may show 
MQCP region when the parameter is tuned away from 
the narrow Ising region. This explains the experimental 
observation on the k— ET compound. 

In this paper, we have considered only MI transitions in 
the ordered phase. We assume that a preexisting gap ex- 
ists in the metallic phase because of the symmetry break- 
ing. However, beyond the mean-field theory, even with- 
out the symmetry breaking, preexisting gap opens in the 
metallic phase as in the case of Mott insulator at low 
dimensions. Therefore, the consequence of the present 
mean-field treatment may survive without a strict order, 
where A may be replaced with the correlation induced 
gap without long-range order. This is indeed corrob- 
orated by various numerical studies performed beyond 
the mean-field approximations, which consistently sup- 
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port the hyperscaling with 5 = 2 and z = 4. 

The present results show that the MI transition is gov- 
erned by the topological change of the Fermi surface, with 
shrinkage (or emergence) at selected momentum points. 
This is different from other types of mean-field theories 
such as the dynamical mean-field approximations, where 
the MI transition is governed instead by the vanishing 
renormalization factor Z. On the verge of the MI tran- 
sition in the metallic side, it possibly means that the 
Fermi surface is reduced to small pockets which violates 
the Luttinger sum rule. If the system undergoes a Lif- 
shitz transition from a large to small Fermi surface in 
the metallic phase, this violation is allowed in the side of 
small pockets. The improvement of the dynamical mean- 
field theory to include the momentum dependence indeed 
suggests the existence of such a Lifshitz transition as we 
discussed in Sec. |U The existence of the small Fermi sur- 
face without folding of the Brillouin zone in the absence 
of the translational symmetry breaking is a fundamen- 
tal issue and should be more thoroughly studied in the 
future. 

The present study has revealed a mechanism of gener- 
ating quantum critical point of the metal-insulator tran- 
sition, which has been unexplored in the literature. The 
marginal quantum critical point emerges because of the 
generation of the gap determined in cooperation with the 
kinetic energy gain of metallic carriers. Its combined and 
nonlinear effect yields pinning of the quasiparticle energy 
in contrast to the rigid band picture. This may be ex- 
pressed by attractive effective interaction of carriers. It 
generates an unconventional universality. 

A remarkable finding is that MQCP is governed by 
the coexistence of quantum fluctuations and the diver- 
gence of the density fluctuations at small wave numbers. 
Such density fluctuations may cause various further and 
deeper effects near MQCP. Its thorough understanding 
is an intriguing issue left for future studies. Among all, 
revealing momentum and frequency resolved density fluc- 
tuations is an extremely important challenging issue of 
experimental studies, which must reveal an anomalous 
and critical enhancement near MQCP as well as near 
the finite-temperature critical point if accurate measure- 
ments are made. Thorough studies of superconducting 
mechanism arising from MQCP [37j are also extremely 
important issue to be explored. 
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APPENDIX A 

In this appendix we show details of the numerical cal- 
culations. In the standard procedure of solving Eqs. (j4"B"|) 
and (|47p , input parameters are effective interaction g and 
particle density n. In this procedure, chemical potential 
fi and the order parameter m are determined by solving 
Eqs. (14*6")) and (|4"T)) self-consistently. This procedure costs 
much computation time, and is not efficient in determin- 
ing [l and m precisely. 

Instead, we use a different procedure. Namely, we take 
A and n as input parameters, while output parameters 
are g and fi. By following this procedure, first, in Eq. 
(|4"T)) for fixed A, we can determine the chemical poten- 
tial \x to satisfy the condition n = 1 independent of Eq. 
([4"6"| . Then, using this fi, we perform the integration of 
the right hand side of Eq. (|46|) and obtain g. Finally, we 
obtain the order parameter m = A/g. Since this proce- 
dure does not require the self-consistent calculations, it 
is possible to calculate Eqs. (|4"6"|) and (I4T[) rapidly and 
precisely. 

Hereafter, we explain the details of integration of Eqs. 
(l46t and (|47|) . To calculate the fc-space integration in 
Eq. (|4T|) . first we fix k y and obtain the particle density 
integrated over k Xl namely n(k y ). It is easy to solve the 
equation 

£ 2 (k X , ky) i ^^{k^kyY + A 2 - (1 = (Al) 

with respect to cos k x for fixed k y analytically since the 
equation is quadratic with respect to cosfe^. Using the 
roots of the Eq. (fAT|) . we obtain n{k y ) exactly. To sim- 
plify the explanation, we only consider the case that the 
two real roots exist for upper band and their absolute 
values are less than 1. In this case, analytically we can 
obtain the location of the point where E + and fi cross, 
namely k r and — fc r [see Fig. |19j. From this, for the up- 
per[lower] band, the particle density for fixed k y is given 
by n{k y )+ = 2k r /2n [n(k y )_ = 2A(fcj,)/27r], where A(A !/ ) 
is the edge of Brillouin zone. The particle density for 
the fixed k y is given as n(k y ) = n(k y ) + + n(k y )-. In 
other cases, in a similar way, we can obtain n(k y ) an- 
alytically. By integrating n(k y ) over k y , we obtain the 
whole particle density n. Along the k y direction, we use 
Newton-Cotes formula and take the number of points in 
the k y direction up to 400000. 

Next, we explain the way of integrating over k in the 
right hand side of Eq. (|4"6"| . We fix k y and perform the 
k x integration first. Along the k x direction, we use dou- 
ble exponential (DE) formula, which is suitable for the 
integration of the analytic function in the integration pe- 
riod. It is known that DE formula is one of the optimized 
methods to integrate the analytic function with the high- 
est accuracy using the smallest number of points. In the 
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above case, integration periods I\ and I2 are given as 
h = [-A, —k r ] and I 2 = [k r ,A] (see Fig. mj). DE for- 
mula is useful in these periods, since the right hand side 
of Eq. (|4l))) is analytical. We take the number of points in 
these periods up to 1000. Then we perform the integra- 
tion over k y using the Newton-Cotes formula, and obtain 
the effective interaction g. We take the number of points 
in the ^-direction up to 400000 for the integration. 
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FIG. 19: (color online) Schematic band structure for fixed k y . 
Circles represent the crossing points of the upper band and 
the chemical potential, which occur at k r . 



where A c = 2t' . Using this, the electron density of the 
upper band, namely electron density X + including spin 
degeneracy is given by 



X 4 



2tt 

1 / 2A C - A + fi 
2^ V A^ 



(B7) 



The location of the top of the lower band is given by 
(7r/2,7r/2) and its equivalent points. Near (7r/2,7r/2), 
using Eq. ((44)) . £L is expanded with respect to q x — 
7r/2 — k x and q y = 7r/2 — k y , which leads to 

E- ~ -A + 2q 2 (t' sin2g - + Sin2g) ) - (B8) 



where we define q = Jq\ + q 2 . Then we obtain 



dE 



ih -4g(— -(*'-—) sin 2^). (B9) 



From Eq. ()B9|) . the DOS of the lower band near 
(7r/2, 7r/2) is given by 



APPENDIX B 

In this appendix we derive Eqs. fM)) . and ([rjt]|) . 
First, we obtain the density of states (DOS) at the bot- 
tom (top) of the upper (lower) band, respectively. 

The location of the bottom of the upper band is given 
by (tt,0) and (0, 7r). Near (ir, 0), using Eq. (gj), E + is 
expanded with respect to q x = k x — ir and q y — k y as 

E + ~ —At' cos q x cos q y 

+ ^4t 2 (cosq x - coaq y ) 2 + A 2 - /i (Bl) 

At' + A + 2t'q 2 - fi, (B2) 



where we define q = yQ x + Q y - From this, we obtain 
dE. 



dq 



Aqt'. 



The DOS of the upper band near (71-, 0) is given by 

1 



Da 



(B3) 



(B4) 



Near the MI transition, the Fermi wave number q + is 
determined by the condition 

E + At' + A + 2t'q\ - (j, = 0, (B5) 



where /i is the chemical potential. From Eqs. (|B2[) and 
(IB5I). we obtain 



2A C - A + fi 



A n 



(B6) 



d6 



(B10) 

Near the MI transition, the Fermi wave number q_ is 
determined by the condition 

£_ ~ -A+2^(t / sin2g- t2(1 + A Sin2g) )~M = 0, (Bll) 



where /Lt is the chemical potential. From Eqs. (|B8p and 
(|Blip . we obtain 



1-{0) 



t'(-n - A) 



t 2 - {2t' z -t 2 )sm26 



(B12) 



Using this, the hole density of the lower band, namely 
the hole density AT_ is given by 



x -= 4 ^ 

(-//-A) 



2tt\A 2 - i' 2 ' 



(B13) 



At half filling, X + should be equal to X^ . From this 
condition, recalling Eq. (fr?3")) and using Eqs. (|B7p and 
(IB13|) . we obtain a as 



2t' - yjt 2 - t' 2 



2t' + yjt 2 - t' 2 



(B14) 
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